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■ Abstract 

We use light-cone QCD sum rules to evaluate the strong coupling g^K+K- 
which enters in several analyses concerning the scalar / (980) meson. The 
result: 6.2 < gf K+K- < 7.8 GeV is larger than in previous determinations. 
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1 Introduction 



Light scalar mesons are the subject of an intense and continuous scrutiny aimed at clari- 
fying several aspects of their nature that still need to be unambiguously established [1, 2]. 
From the experimental point of view, these particles are difficult to resolve because of the 
strong overlap with the continuum background. On the other hand, the identification is 
made problematic since both quark-antiquark (qq) and non qq scalar states are expected 
to exist in the energy regime below 2 GeV. For example, lattice QCD and QCD sum 
rule analyses indicate that the lowest lying glueball is a ++ state with mass in the range 
1.5-1.7 GeV [3]. Actually, the observed light scalar states are too numerous to be acco- 
modated in a single qq multiplet, and therefore it has been suggested that some of them 
escape the quark model interpretation. In addition to glueballs, other interpretations 
include multiquark states and admixtures of quarks and gluons. 

Particularly debated is the nature of the meson / (980). Among the oldest suggestions, 
there is the proposal that quark confinement could be explained through the existence 
of a state with vacuum quantum numbers and mass close to the proton mass [4]. On 
the other hand, following the quark model and considering the strong coupling to kaons, 
/o(980) could be interpreted as an ss state [5, 6, 7, 8]. However, this does not explain 
the mass degeneracy between / (980) and ao(980) interpreted as a (uu — dd) / \/2 state. A 
four quark qqqq state interpretation has also been proposed [9]. In this case, /o(980) could 
either be nucleon-like [10], i.e. a bound state of quarks with symbolic quark structure 
fo = ss(uu + dd)/\/2, the ao(980) being ao = ss(uu — dd)/y/2, or deuteron-like, i.e. a 
bound state of hadrons. If fo is a bound state of hadrons, it is usually referred to as a KK 
molecule [11, 12, 13, 14]. In the former of these two possibilities the mesons are treated 
as point-like, while in the latter they should be considered as extended objects. 

The identification of the fo and of the other lightest scalar mesons with the Higgs nonet 
of a hidden U(3) symmetry has also been suggested [15]. Finally, a different interpretation 
consists in considering /o(980) as the result of a process in which strong interaction en- 
riches a pure qq state with other components, such as KK^j, a process known as hadronic 
dressing [6, 16]; such an interpretation is supported in [2, 5, 6, 8, 17, 18, 19]. In ref. [20] 
it has been shown that the experimentally observed lightest scalar particles in the 1=1 
and 1=1/2 sectors can be reproduced in this way, starting from a bare qq and sq structure 
respectively (q being a light non strange quark). On the other hand, 1=0 states are the 
most elusive ones, since there are two possible bare structures, qq and ss, which could 
not only undergo hadronic dressing, but also mix through hadronic loops. The resulting 
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picture strongly depends on the couplings of the bare structures to the hadronic channels. 

Several experimental analyses aimed at discriminating among the different possibili- 
ties. In particular, the radiative — > f Q j decay mode has been identified as an effective 
tool for such a purpose [10, 12, 21]. As a matter of fact, if f has a pure strangeness 
component fo = ss, the dominant — > /o7 decay mechanism is the direct transition, 
while in the four-quark scenario <fi — > /o7 is expected to proceed through kaon loops with 
a branching fraction depending on the specific bound state structure [12, 21]. 

An important hadronic parameter entering in several analyses involving / (980) is the 
strong coupling gf K+K-- Indeed, the kaon loop diagrams contributing to <fi — > / 7 are 
expressed in terms of gf K+K~, as well as in terms of the coupling g^K+K- which can 
be inferred from experimental data on meson decays. The coupling gf K+K- can be 
obtained from various processes, and we shall present an overview of the determinations 
in the last part of this paper. It is interesting to carry out a calculation in a framework 
based on QCD, trying to point out what is a distinctive feature of the scalar particles, 
i.e. their large couplings to the hadronic states. 

The present study is devoted to a determination of gj K+K- by light-cone QCD sum 
rules, a method applied to the calculation of several hadronic parameters both in the 
light, both in the heavy quark sector [22, 23]. The analysis and the numerical results are 
presented in Section 2, while a summary of the experimental data and of other theoretical 
determinations is given in Section 3. 

2 Coupling gf K+K- by light-cone QCD sum rules 

In order to evaluate the strong coupling gf K+K~, defined by the matrix element: 

(K + (q)K-(p)\f (p + q))=g foK+K ~ , (1) 
we consider the correlation function 

T,(p,q) = i J d 4 xe^ (K + (q)\T[J«(x)J fo (0W) • (2) 

The quark currents and J/ represent the axial- vector = ^7^75 s and the scalar 
Jf = ss current, respectively, while the external kaon state has four momentum q, with 
q 2 = M\. The choice of the J/ = ss current does not imply that / (980) has a pure ss 
structure, but it simply amounts to assume that J/ has a non-vanishing matrix element 
between the vacuum and f [19, 24]. Such a matrix element, as mentioned below, has 
been derived by the same sum rule method. 
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Exploiting Lorentz invariance, T^ can be written in terms of two independent invariant 
functions, Ti and T 2 : 

T>, q) = i T.ip 2 , (p + qf) p, + T 2 (p 2 , (p + g) 2 ) g M . (3) 

The analysis of the correlation function in eq.(2), following the general strategy of QCD 
sum rules, allows us to obtain a quantitative estimate of gf K+K-- The method consists 
in representing T M in terms of the contributions of hadrons (one-particle states and the 
continuum) having non-vanishing matrix elements with the vacuum and the currents 
Jff and Jy , and matching such a representation with a QCD expression computed in a 
suitable region of the external momenta p and p + q [25]. 

Let us consider, in particular, the invariant function Ti in eq.(3) that can be repre- 
sented by a dispersive formula in the two variables p 2 and (p + q) 2 : 

The hadronic spectral density p had gets contribution from the single-particle states K and 
/o, for which we define current-particle matrix elements: 

(/o(980)(p + 9 )|J /o |0) = M /o / (5) 
(0\J?\K{p))=if K p ll , (6) 

as well as from higher resonances and a continuum of states that we assume to contribute 
in a domain D of the s, s' plane, starting from two thresholds Sq and s' . Therefore, 
neglecting the fo width, the spectral function p had can be modeled as: 

P had (s, s') = f K M f Jg foK+K -5(s - M 2 K )5(s> - M 2 ,) + p cont (s, s>)6(s - s )6(s> - s' Q ) , (7) 

where p cont includes the contribution of the higher resonances and of the hadronic contin- 
uum. The resulting expression for 7\ is: 

T, (V 2 (V + q) 2 ) = fK M fof9fo* + K- + f dsds > P C ° n \s,s') 

1 {P ' {P + q) } (Mi - p=») {M % -{p + qf) + Id (s - p*) [s' -(p + q) 2 ]- ( " j 

We do not consider possible subtraction terms in eq.(4) as they will be removed by a 
Borel transformation. 

For space-like and large external momenta (large — p 2 , —{p + q) 2 ) the function Ti 
can be computed in QCD as an expansion near the light-cone x 2 = 0. The expansion 
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involves matrix elements of non-local quark-gluon operators, which are defined in terms of 
kaon distribution amplitudes of increasing twist. 1 The first few terms in the expansion 
are retained, since the higher twist contributions are suppressed by powers of l/(— p 2 ) 
or l/(— (p + q) 2 ). As a result, the following expression for 7\ is obtained to twist four 
accuracy: 

1 
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The functions ip p and ip a appearing in eq.(9) are kaon distribution amplitudes defined by 
the matrix elements 

f M 2 /■! 

< K(q)\u(x)i l5 s(0)\0 > = du e iuq - x <p p (u) , (10) 

TDj g Jo 

<K(q)\u(x)a^ 5 s(0)\0> = i{q^x v - q vX J»L-JL J du e^ x Va (u) (11) 



m s being the strange quark mass (we put to zero the mass of the light quarks). 2 Moreover, 
g2{u) is defined by the matrix element 

< K(q)\u(x) lfll5 s(0)\0 >= - if K% f 1 du e iu * x [y K {u) + x 2 9l (u)} 

J 

+ f K (x, - I du e iu ^g 2 (u) . (12) 

q ■ x Jo 

Kaon matrix elements of quark-gluon operators also contribute to eq.(9); they are param- 
eterized in terms of twist three and twist four distribution amplitudes: 

< K(q)\u(x)a a p-f 5 g s G^(vx)s(0)\0 >= 

if3K[(q^q a g^ - qvq a g^) - (q^g^a - q^g^a)] J £>o; (p 3K (ai)e iq ' x{ai+va3 \ (13) 



1 The short-distance expansion of the 3-point vacuum correlation function of one scalar ss and two 
pseudoscalar sij^q densities has been considered in [26]. The present calculation mainly differs for the 
possibility of incorporating an infinite series of local operators [23] . 

2 Thc path-ordered gauge factor Pexpig s J* dtx 11 A^{tx) is not included in the matrix elements having 
chosen the gauge x^A^ = 0. 
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The operator G a/ 3 is the dual of G a p\ G a p = \e a p&pG 5p ; Va>i is defined as Voti = 
da\da2da 3 5{l — a±— 0:2 — 0:3). The function (p 3 x is twist three, while the distribution ampli- 
tudes in (14) and (15) are twist four. The functions tp and appearing in eq.(9) are defined 
in terms of ip±, <p\\, (p±_ and <p\\ as follows: ^{a 3 ) = — Jq 3 dt J 1_< da^(ai, 1 — a± — t,t), 
and 0(o;j) = - Jq 1 dt$(t, 1 - t - a 3 , a 3 ), with $ = ip± + (p\\ + <p± + <p\\. 

The sum rule for gf K+K- follows from the approximate equality of eqs.(8) and (9). 
Invoking global quark-hadron duality, the contribution of the continuum in (8) can be 
identified with the QCD contribution above the thresholds Sq, s' . This allows us to iso- 
late the pole contribution in which the coupling appears. The matching between the 
expressions in (8) and (9) can be improved performing two independent Borel transfor- 
mations with respect to the variables — p 2 and — (p + q) 2 . Defining M\ and M| as the 
Borel parameters associated to the channels p 2 and (p + q) 2 , respectively, and using the 
identity: 
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, we get the following expression for the Borel 
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where ml = m 2 + u (l — u )M^. Analogously, a double Borel transformation can be 
carried out for the hadronic representation eq.(8): 
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As shown by (18), the Borel transformation exponentially suppresses the contribution 
of the higher states and of the continuum; furthermore, it removes possible subtraction 
terms in (4) which depend only on p 2 or {p + q) 2 . 

The second term in (18) represents the continuum contribution. In order to identify 
it with part the QCD term (17), a prescription has been proposed in [27]. It consists 
in considering the symmetric points M 2 = Mf = 2M 2 (corresponding to m = 1/2) in 
the (M 2 ,Mf) plane and performing the continuum subtraction through the substitution 

.2 - 2 

e ~m* ^ e ~m* — eM^in the leading-twist term in (17). Such a prescription is not 
adeguate in our case, where the Borel parameters correspond to channels with different 
mass scales and should not be constrained to be equal. Here we can exploit the property 
of the amplitudes <p p (u) and (p a (u) of being polynomials in u (or 1 — u): 

1 N 

Vp (u) + - V ' a {u) = Y J b k ^-u) k (19) 

fc=0 

in order to compute their contribution in the duality region D. As for the other terms in 
(17), they represent a small contribution to the QCD side of the sum rule, and therefore 
the calculation can leave them unaffected. 
The final expression for gf K+K- reads: 
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with A = S ° M " S and s the smallest continuum threshold. The prescription in [27] is 
obtained for M 2 = Mf , % — and neglecting terms of order M\. An interesting feature 



of eq.(20) is that, changing M 2 and Mf independently, it is possible to vary the point 
u where the distribution amplitudes are evaluated and contribute, while in the standard 
approach the final result is essentially related to the value of the distribution amplitudes 
in a selected point. 

The main nonperturbative quantities constituting the input information in the sum 
rule (20) are the kaon light-cone wave functions. A theoretical framework for their deter- 
mination relies on an expansion in terms of matrix elements of conformal operators [28]. 
For the function ip p , conformal expansion results in the expression 

<p p (u,n) = £«SMCj?(0 (21) 
k 

with £ = 2u — 1, Oq = 1 and C[ the Gegenbauer polynomials. In (21) we have included 
the normalization scale dependence of the distribution amplitude <p p , which appears in 
the multiplicatively renormalizable coefficients a p k {p). The nonperturbative information is 
encoded in the coefficients, which are peculiar for the various mesons. In the case of kaon, 
the asymmetry between the strange and nonstrange quark momentum distribution in the 
meson can be taken into account by non- vanishing odd-order coefficients a v k . Such SU(3) 
flavour violating effects have not been investigated so far for distribution amplitudes of 
twist larger than two, and we neglect them in the following, with consequences that we 
shall mention below. As for the even order coeficients, their updated values are reported 
in [27, 29]: a\ = 30r] 3 - \p 2 and a\ = -3rj 3 u 3 - §p 2 - %p 2 a 2 , with a 2 = 0.2, r] 3 = 0.015, 
u 3 = —3 at the scale // ~ 1 GeV. We have taken into account the meson mass corrections, 

2 

related to the parameter p 2 = ^k-, worked out in [29]. 

K 

Analogously, the ip a distribution amplitude can be expressed as 

<Pa(u,n) = Ml -«)E a *M<?*(0 ( 22 ) 

k 

with dp = 1, a 2 = 5?73 — ^i] 3 cu 3 — ^p 2 — ^p 2 a 2 . For p 3 K(cti) and for the other higher twist 
distribution amplitudes we refer to the expressions reported in [27, 29]. 

In the analysis of eq.(20) we use m s {lGeV) = 0.14 GeV [30], M K = 0.4937 GeV, 
M fo = 0.980 GeV, f K = 0.160 GeV and / = (0.180 ± 0.015) GeV [19]. The threshold 
parameter s is varied around the value s = 1.1 GeV 2 fixed from the determination of 
fx using two-point sum rules [31]. 

The result for gf x+K- versus the Borel parameters M 2 and Mf is depicted in fig.l. A 
stability region where the outcome does not depend on M 2 can be selected. Such a region 
does not correspond to the line Mf = Mf , but to the range 0.8 < M 2 < 1.6 GeV 2 with 
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Figure 1: Coupling gf K+K- as a function of the Borel parameters Mf and Mf , for so = 1.1 
GeV 2 . 

Mf extending up to Mf ~ 5 GeV 2 . Varying Mf and Mf in this region, and changing 
the values of the thresholds and of the other parameters, we obtain the result depicted in 
fig.2, which can be quoted as 6.2 < gf K+K- < 7.8 GeV. 

Let us briefly discuss the uncertainties affecting the numerical result. As for the 
SU(3)p breaking effects rendering the kaon distribution amplitudes asymmetric with re- 
spect to the middle point, the neglect should have a minor role in our approach, due to 
the possibility of exploring wide ranges of the variable u and smoothing the effects of the 
actual shapes of the wave functions. Another uncertainty is related to the value of the 
strange quark mass, m s ; since the dependence of the sum rule on m s mainly involves the 
ratio M\jm ai one can fix this ratio using chiral perturbation theory, obtaining results in 
the same range quoted for gj K+K-- 

We can compare now our result with the available experimental determinations of 
9f K+K-, as we h as with the results of other calculations. We shall see how complex the 
scenario is. 




Figure 2: Coupling gf x+K- as a function of the Borel parameter Mf , varying s in the 
range 1.05 < s < 1.15 GeV 2 and M 2 in the range 0.7 < M\ < 2.0 GeV 2 . 



3 Other determinations of gf 



K+K- 

As discussed in the Introduction, g/ K+K- can be considered in connection with the ra- 
diative (ft — > /07 decay mode. Several analyses go through this decay channel. KLOE Col- 
laboration at the DA$NE collider in Frascati has examined the decay channel — > 7r 7T°7 
measuring the branching fraction: B(<f> — > 7r°7r°7) = (1.09 ± 0.03 stat ± 0.05 syst ) x 10 -4 
[32]. The decay mode is supposed to proceed through pit intermediate state and through 
kaon loop processes, with the kaons annihilating into scalar resonances that subsequently 

decay to 7r°7r°. Different fits of the two pion invariant mass spectrum are per- 

formed in order to measure the parameters of the scalar states. In a first fit (A) only 
the contribution of the intermediate state / (980) is considered, and the three parameters 
,k+k~ an d 9 2 f K+K-/9 2 f 07 nr are determined. In a second fit (B) the contribution 
of a possible broad scalar a state is included, and the coupling is considered as a 
further parameter. It is assumed that the two pion decay modes saturate the / width, 
and that B(fo — * 7i + 7i~) = 2 B(fo — > 7r°7r°) invoking isospin symmetry. Fit A provides 

B(<f> -> /o7) = (3.3 ±0.2) x 10- 4 and 9f ° K+I <~ = 1.29 ±0.14 GeV 2 (x 2 /ndf = 109.53/34). 

(4tt) 
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Fit B gives instead: BU -> / t) = (4.47 ±0.21) x 10~ 4 and 9f f +K = 2.79 ±0.12 GeV 2 

(4tt) 

(x 2 /ndf = 43.15/33). The negative interference between the contributions of the broad 
a and the /o is responsible of the improvement in the accuracy of the fit. In both cases 
sizeable values for gf x+K- are obtained; they are reported in Table 1. 

An analogous analysis has been performed by the CMD-2 Collaboration at the VEPP- 
2M collider in Novosibirsk. From a combined fit to the spectra of the decays — > 7r + 7r - 7 
and (f) -> 7r°7r°7, CMD-2 Collaboration obtains: B{(j> -> /07) = (2.90 ±0.21 ± 0.65) x 10~ 4 

and ^7^p = 1.48±0.32 GeV 2 [33]. A similar result is quoted by the SND Collaboration 

Q 2 

at the same VEPP collider: B(<f> -> /07) = (3.5 ± 0.3d=J;|) x 10~ 4 and = 
2A7toil GeV 2 [34]. 

Other determinations of gf K+K- rer y 011 the analysis of different physical processes. 
Considering the central fo production in collisions, the WA102 experiment at CERN 

g 2 + - 

gets: |^ j — = 0.38 ± 0.06 GeV 2 [35]. On the other hand, analyzing the f production 
in D s decays to three pions, the Collaboration E791 at Fermilab finds a value compatible 
with zero [36]. These results are also reported in Table 1. 

In Ref.[10] the coupling constant is evaluated for different values of the phase shift of 
the elastic background in the ttty — > nn reaction, of the ratio R = g 2 [)K + K - /g 2 o7T+7T - and 
according to different scenarios for the fo structure, obtaining results in a wide range: 
9f oK+ K- e [1.95,7.3] GeV. 

The analysis of the decay channel J/ip — > 4>KK(ix-k) has been carried out in Ref.[37]. 
The fo pole is described as a Breit-Wigner resonance coupled to two channels. Two fits 
of the experimental data are performed depending upon the mi phase shift data used, 
obtaining g j k+k- — 2-5 ± 0.15 GeV and g^K+K- — 2-0 ± 0.06 GeV, respectively. 

A prediction for gf K+K- based on chiral symmetry and the linear sigma model, when 
no mixing with the a is considered, is: g/ K+K- = 2.24 GeV [38], to be compared to old 
determinations gj K+K- = 2.74 GeV [39]. Using the method of the T-matrices, the value 
9f K+K- = 3.8 GeV is obtained [40]. 

Considering all the above results one sees that a general consensus on gf x+K- has not 
been reached, so far. In particular, experimental analyses of different processes produce 
contradicting results. The outcome from <fi — > / 7 points towards sizeable values of the 
coupling, consistent with the light-cone sum rule result. One has to say that the error 
quoted for the experimental determinations, which in general looks small, mainly accounts 
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for the statistical uncertainties; one could infer the size of the systematical uncertainties 
comparing different determinations. 

As for gf K+K- from D s decays, presumably the determination will be improved at 
the B factories by experiments such as BaBar at SLAC, thanks to large available samples 
of D s mesons. In these measurements gf K+K- is expected to be determined by coupled 
channel analyses, with D s decaying to final states containing kaons as well as pions [41]. 

Table 1: Experimental determinations of g/ K+K- using different physical processes. Dou- 
ble items refer to two different fits (see text). 



Collaboration 


process 


9f K+K- (GeV) 


Ref. 


KLOE 


-> fol (A) 
<P^fol(B) 


4.0 ± 0.2 (A) 
5.9 ± 0.1 (V) 


[32] 


CMD-2 


4> -»• hi 


4.3 ±0.5 


[33] 


SND 


-»• foi 


5.6 ±0.8 


[34] 


WA102 


pp 


2.2 ±0.2 


[35] 


E791 


D s -> 3tt 


0.5 ±0.6 


[36] 



4 Conclusions 

The purpose of this paper was the evaluation of the strong coupling constant g/ K+K-, the 
value of which is rather controversial, as it emerges comparing different experimental and 
theoretical determinations. In particular, the KLOE Collaboration measured a larger 
value than in other determinations, with a greater accuracy as well. However, such a 
result stems from the investigation of <fi — > fol, and therefore it is mandatory to wait 
for the study of unrelated processes, namely the combined analysis of D s decays to pions 
and kaons. The outcome of light-cone QCD sum rules is in keeping with a large value for 
the coupling. The uncertainty affecting the result is intrinsic of the method and does not 
allow a better comparison with data. However, the analysis confirms a peculiar aspect of 
the scalar states, i.e. their large hadronic couplings, thus pointing towards a scenario in 
which the process of hadronic dressing is favoured. 

Acknowledgments We thank A. Palano, T. Aliev and M.E. Boglione for discussions. 
We acknowledge partial support from the EC Contract No. HPRN-CT-2002-00311 (EU- 
RIDICE). 



12 



References 

[1] For reviews see: L. Montanet, Rep. Prog. Phys. 46 (1983) 337; F.E. Close, Rep. 
Prog. Phys. 51 (1988) 833; N.N. Achasov, Nucl. Phys. Proc. Suppl. 21 (1991) 189; 
M. R. Pennington, Proceedings of HADRON '95, M.C. Birse, G.D. Lafferty, J.A. 
McGovern eds., World Scientific (1996) page 3; T. Barnes, hep-ph/0001326. 

[2] F. E. Close and N. A. Tornqvist, J. Phys. G 28 (2002) R249. 

[3] C. J. Morningstar and M. J. Peardon, Phys. Rev. D 60 (1999) 034509; S. Narison, 
Nucl. Phys. A 675 (2000) 54C. 

[4] F. E. Close, Y. L. Dokshitzer, V. N. Gribov, V. A. Khoze and M. G. Ryskin, Phys. 
Lett. B 319 (1993) 291. 

[5] N.A. Tornqvist, Phys. Rev. Lett. 49 (1982) 624. 

[6] N.A. Tornqvist, Z. Phys. C 68 (1995) 647. 

[7] N.A. Tornqvist and M. Roos, Phys. Rev. Lett. 76 (1996) 1575. 

[8] E. van Beveren et al, Z. Phys. C 30 (1986) 615; M.D. Scadron, Phys. Rev. D 26 
(1982) 239; E. van Beveren, G. Rupp and M.D. Scadron, Phys. Lett. B 495 (2000) 
300. 

[9] R.L. Jaffe, Phys. Rev. D 15 (1977) 267, 281; D 17 (1978) 1444; R.L. Jaffe and K. 
Johnson, Phys. Lett. B 60 (1976) 201. 

[10] N.N. Achasov and V.N. Ivanchenko, Nucl. Phys. B 315 (1989) 465; N.N. Achasov 
and V.V. Gubin, Phys. Rev. D 56 (1997) 4084. 

[11] J. Weinstein and N. Isgur, Phys. Rev. Lett. 48 (1982) 659; Phys. Rev. D 27 (1983) 
588; Phys. Rev. D 41 (1990) 2236. 

[12] N. Brown and F.E. Close, in The DA&NE Physics HandBook, L. Maiani, G. Pancheri 
and N. Paver eds, INFN Frascati, 1995, pp. 447-464; F.E. Close, N. Isgur and S. 
Kumano, Nucl. Phys. B 389 (1993) 513. 

[13] R. Kaminski, L. Lesniak and J.P Maillet, Phys. Rev. D 50 (1994) 3145. 

[14] N.N. Achasov, V.V. Gubin and V.I. Shevchenko, Phys. Rev. D 56 (1997) 203. 



13 



[15] N. A. Tornqvist, arXiv:hep-ph/0204215. 

[16] M. Boglione and M. R. Pennington, Phys. Rev. Lett. 79 (1997) 1998. 

[17] N. A. Tornqvist, arXiv:hep-ph/0008136. 

[18] C. M. Shakin and H. Wang, Phys. Rev. D 63 (2001) 014019. 

[19] F. De Fazio and M. R. Pennington, Phys. Lett. B 521 (2001) 15. 

[20] M. Boglione and M. R. Pennington, Phys. Rev. D 65 (2002) 114010. 

[21] F. E. Close and A. Kirk, Phys. Lett. B 515 (2001) 13. 

[22] I. I. Balitsky, V. M. Braun and A. V. Kolesnichenko, Nucl. Phys. B 312 (1989) 509; 
V. L. Chernyak and I. R. Zhitnitsky, Nucl. Phys. B 345 (1990) 137. 

[23] For a recent review of the method see: P. Colangelo and A. Khodjamirian, in "At the 
Frontier of Particle Physics", vol. 3, M. Shifman ed., World Scientific (2001), page 
1495. 

[24] T. M. Aliev, A. Ozpineci and M. Savci, Phys. Lett. B 527 (2002) 193. 

[25] T-products of quark currents between the vacuum and one-particle states were first 
investigated in N. S. Craigie and J. Stern, Nucl. Phys. B 216 (1983) 209. 

[26] S. Narison, Nucl. Phys. B 509 (1998) 312. 

[27] V. M. Belyaev, V. M. Braun, A. Khodjamirian and R. Riickl, Phys. Rev. D 51 (1995) 
6177. 

[28] V. M. Braun and I. E. Halperin, Z. Phys. C 48 (1990) 239. 
[29] P. Ball, JHEP 9901 (1999) 010. 

[30] P. Colangelo, F. De Fazio, G. Nardulli and N. Paver, Phys. Lett. B 408 (1997) 340; 
K. Maltman and J. Kambor, Phys. Rev. D 65 (2002) 074013; M. Jamin, J. A. Oiler 
and A. Pich, Eur. Phys. J. C 24 (2002) 237. 

[31] A. A. Ovchinnikov and A.Pivovarov, Phys. Lett. B 163 (1985) 231. 

[32] A. Aloisio et al. [KLOE Collaboration], Phys. Lett. B 537 (2002) 21. 

14 



[33] R. R. Akhmetshin et al. [CMD-2 Collaboration], Phys. Lett. B 462 (1999) 380; Nucl. 
Phys. A 675 (2000) 424C. 

[34] M. N. Achasov et al, Phys. Lett. B 485 (2000) 349. 

[35] D. Barberis et al. [WA102 Collaboration], Phys. Lett. B 462 (1999) 462. 

[36] E. M. Aitala et al. [E791 Collaboration], Phys. Rev. Lett. 86 (2001) 765. 

[37] R. Escribano, A. Gallegos, J. L. Lucio M, G. Moreno and J. Pestieau, arXiv:hep- 
ph/0204338. 

[38] M. Napsuciale, arXiv:hep-ph/9803396; J. L. Lucio Martinez and M. Napsuciale, Phys. 
Lett. B 454 (1999) 365. 

[39] J. F. Donoghue, B. R. Holstein and G. Valencia, Phys. Rev. D 35 (1987) 2769; 
G. Ecker, A. Pich and E. de Rafael, Phys. Lett. B 189 (1987) 363; G. Ecker, A. Pich 
and E. de Rafael, Nucl. Phys. B 303 (1988) 665; L. M. Sehgal, Phys. Rev. D 38 (1988) 
808; S. Nussinov and T. N. Truong, Phys. Rev. Lett. 63 (1989) 1349 [Erratum-ibid. 
63 (1989) 2002]; J. L. Lucio Martinez and J. Pestieau, Phys. Rev. D 42 (1990) 3253. 

[40] J. A. Oiler and E. Oset, Nucl. Phys. A 620 (1997) 438 [Erratum-ibid. A 652 (1999) 
407]; J. A. Oiler, arXiv:hep-ph/0205121. 

[41] B. T. Meadows, arXiv:hep-ex/0210065. 



15 



